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We present a detailed study of chemical freeze-out in nucleus-nucleus collisions at beam ener- 
gies of 11.6, 30, 40, 80 and 158^1 GeV. By analyzing hadronic multiplicities within the statistical 
hadronization approach, we have studied the strangeness production as a function of centre of mass 
energy and of the parameters of the source. We have tested and compared different versions of 
the statistical model, with special emphasis on possible explanations of the observed strangeness 
hadronic phase space under-saturation. We show that, in this energy range, the use of hadron 
yields at midrapidity instead of in full phase space artificially enhances strangeness production and 
could lead to incorrect conclusions as far as the occurrence of full chemical equilibrium is concerned. 
In addition to the basic model with an extra strange quark non- equilibrium parameter, we have 
tested three more schemes; a two-component model superimposing hadrons coming out of single 
nucleon-nucleon interactions to those emerging from large fireballs at equilibrium, a model with 
local strangeness neutrality and a model with strange and light quark non-equilibrium parameters. 
The behaviour of the source parameters as a function of colliding system and collision energy is 
studied. The description of strangeness production entails a non-monotonic energy dependence of 
strangeness saturation parameter 7s with a maximum around 30A GeV. We also present predic- 
tions of the production rates of still unmeasured hadrons including the newly discovered 0^(1540) 
pentaquark baryon. 



The main goal of the ultra-relativistic nucleus-nucleus (A- A) collisions programme is to create in terrestrial labora- 
tories a new state of matter, the Quark-Gluon Plasma (QGP). The existence of this phase, where quarks and gluons 
are deconfined, i.e. can freely move over several hadronic distances, is a definite prediction of quantum chromodynam- 
ics (QCD). In a search for QGP signals A-A collisions at different centre of mass energies per nucleon-nucleon (NN) 
pair have been studied: from few GeV to several hundreds of GeV recently attained in Au-Au collisions at RHIC. 

Recently, accurate measurements of hadron production in central Pb-Pb collisions at 40, 80 and 158^ GeV of beam 
energy became available U and also preliminary data at iOA GeV have been presented following an energy scan 
programme carried out by the experiment NA49 at CERN SPS. This programme is motivated by the hypothesis Q 
that the threshold for creation of QGP in the early stage of Pb-Pb collisions might be located in the low SPS energy 
range, roughly between 20 and 4QA GeV of beam energy. 

One of the main results of the study of high energy A-A collisions is a surprising success of the statistical-thermal 
models in reproducing essential features of particle production 0, 0, IE S tS,, A 10.. .11.1 . This model succeeds also 
in describing particle multiplicities in many kinds of elementary collisions [12, HJIi suggesting that statistical 
production is a general property of the hadronization process itself |l3lll5j |. Furthermore, the statistical hadronization 
model (SHM) supplemented with the hydrodynamical expansion of the matter, to a large extent also reproduces 
transverse momentum spectra of different particle species [iS]- 

Hence, the SHM model proves to be a useful tool for the analysis of soft hadron production and particularly to 
study strangeness production, whose enhancement has since long been proposed as a signature of QGP formation. 



I. INTRODUCTION 
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Furthermore, anomalies in the energy dependence of strangeness production have been predicted as a signature of 
deconfinement and have been indeed observed experimentally jl7j , suggesting that the onset of the phase transition 
could be located around 30A GeV. It is thus important to make a systematic analysis, within the framework of SHM, 
of the presently available hadronic multiplicities measured in Pb-Pb collisions at 30, 40 and 80A GeV, which - to our 
knowledge - is done here for the first time. 

Along with these intriguing questions, our work is also motivated by issues related to the application of statistical 
model itself. In fact, different versions of this model have been used in the past by different authors leading to 
somewhat different results and conclusions. These mainly stem from the alternative use of midrapidity and full phase 
space multiplicities, from the allowance of non-equilibrium abundances of hadrons, from the assumption of exact local 
vanishing strangeness etc. Therefore, we consider the comparison of these different approaches a worthwhile step. 
This has been made it possible by now by the availability of an accurate and large multiplicity sample in Pb-Pb 
collisions at 158A GeV as well as the corresponding data for pp interactions. 

The paper is organized as follows: a brief description of our main version of the SHM is given in Sect. 2. The 
experimental data selected for the analysis are summarized in Sect. 3. In Sect. 4 the results of the analysis using the 
main version and alternative schemes of the SHM are given. Finally, in Sect. 5 we present and discuss the energy 
dependence of the chemical freeze-out stage. Summary and conclusions are drawn in Sect. 6. 

II. THE STATISTICAL HADRONIZATION MODEL 

The main idea of the SHM is that hadrons are emitted from regions at statistical equilibrium, called clusters or 
fireballs. No hypothesis is made about how statistical equilibrium is achieved; this can be a direct consequence of the 
hadronization process. In a single collision event, there might be several clusters with different collective momenta, 
different overall charges and volumes. However, Lorentz- invariant quantities like particle multiplicities are independent 
of clusters momenta, while they depend on charges and volumes. If final state interactions among formed hadrons 
occur, particle multiplicities are frozen when inelastic interactions cease (chemical freeze-out). Thus, by analyzing 
measured hadron abundances, a snapshot is taken of clusters at that particular stage of the evolution, which may 
significantly precede the final kinetic freeze-out stage, when also elastic interactions cease. However, it should be 
pointed out that chemical and kinetic freeze-out may depend on the hadron species and the assumption of a single 
chemical freeze-out is certainly an approximation. Most calculations in SHM are carried out in the framework of the 
ideal hadron-resonance gas, that is handling resonances as free particles: this amounts to take a considerable part of 
the hadronic interactions between strongly stable hadrons into account Ts^ . 

As has been mentioned, final multiplicities depend on the distribution of initial conserved charges (baryon number, 
strangeness and electric charge) among the produced clusters. This distribution is determined by the dynamics of the 
collision and is thus needed as an external input to the statistical model. However, most analyses, including ours, are 
carried out by assuming a single fireball. This is possible provided that one of the two conditions below is fulfilled: 

1. all clusters are large enough to allow a grand-canonical description and all of them have the same values of 
relevant intensive parameters, i.e. temperature and chemical potentials; 

2. clusters are small and must be treated canonically (i.e. counting those states having exactly the same charges 
as the cluster itself) , yet they have the same temperature and the distribution governing fluctuations of charges 
is the same as that obtained by splitting one large cluster - the equivalent global cluster EGG - having as volume 
the sum of all clusters rest frame volumes and charges the sum of all clusters charges (see Appendix A). In this 
case the overall particle multiplicities turn out to be those calculated in the canonical, perhaps grand-canonical, 
ensemble of the equivalent global cluster • The reduction to EGG could be achieved even for micro-canonical 
clusters with additional requirement on mass fluctuations [T^ . 

The first condition sets stronger requirements and applies in the Bjorken's boost-invariant scenario, where all clusters 
are to have the same parameters independently of their rapidity. The second condition is altogether weaker and 
leaves room for the compatibility between the single fireball analysis and a variation of net baryon number density 
in rapidity. This has been discussed in detail in ref. 0- The argument can be summarized as follows: particle 
multiplicities, being Lorentz invariants, are unaffected by a shift in rapidity of the clusters; therefore, clusters arising 
from the splitting of the EGG can be ordered in rapidity according to their net baryon number without affecting 
fully integrated particle multiplicities and, at the same time, giving rise to an effective variation of the baryon density 
profile. Although the second condition is certainly more appropriate in the examined energy range, from AGS to 
SPS, it must be pointed out that this should not be expected to precisely match physical reality, as well as the first 
condition in its domain of applicability. In other words, discrepancies (hopefully small) between calculations based 
on this model and measurement are to be expected, so that these analyses shall not provide perfect fits even though 
the statistical model was the underlying true model. 
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In this paper we will stick to the picture outlined in the second condition, which implicitely requires the use of full 
phase space multiplicities in order to (hopefully) integrating out correlations between clusters' momenta and charges. 
Besides their general fitness, iir multiplicities also allow to safely enforce overall strangeness neutrality. As has been 
mentioned, if the second condition applies, the multiplicity of any hadron j can be calculated in the canonical ensemble 
of the EGC. Hence, as the EGG has a much larger volume than single clusters', the grand-canonical ensemble, where 
charges are conserved on average, can be a good approximation (see Appendix A). This is the case for the collisions 
examined in this paper [19|. In this case the mean primary multiplicity of the j'^ hadron with mass rrij and spin Jj 
reads: 

where T is the temperature, V the EGG volume, — {Qj,Bj,Sj) is a vector having as components the electric 
charge, baryon number and strangeness of the hadron and /x — {fJ-q, fis, ^J■s) is a vector of the corresponding chemical 
potentials; the upper sign applies to fermions, the lower to bosons. In order to correctly reproduce the data, it is 
also necessary to introduce at least one non-equilibrium parameter suppressing hadrons containing valence strange 
quarks, 7s 7^ 1 |2fll |. With this supplementary parameter, hadron multiplicity is as in Eq. with the replacement: 

exp[/x • qj/T] cxp[/x • qj/T]-/s'' (2) 

where ris stands for the number of valence strange quarks and anti-quarks in the hadron j. 

The abundances of resonances is calculated convoluting |^ with a relativistic Breit-Wigner distribution over a mass 
interval [m — 5m, m + 5m], where 5m — min[m — mthroshoid, 2r]. The minimum mass mthrcshoid is needed to open 
all decay modes. Finally, the overall multiplicity to be compared with the data, is calculated as the sum of primary 
multiplicity ^ and the contribution from the decay of heavier hadrons: 

{n,) = (n,;)P~ + J2 - j){nk), (3) 
k 

where the branching ratios are taken from the latest issue of the Review of Particle Physics '2l'| and the summation 
runs over decays which contribute to the experimentaly measured multiplicity. Among the hadrons and resonances 
contributing to the sum in Eq. ||2Jl, in this work all known states quoted in ref. up to a mass of 1.8 GeV are 
included (see discussion in Sect. 4). 

What we have hitherto described is the main version of the SHM used for the data analysis, that will be henceforth 
referred to as SIIM(75). As has been mentioned in the Introduction, in this work we also test other schemes and 
versions of the SHM, which will be described in detail in Sect. 4. 

III. EXPERIMENTAL DATA SET 

The bulk of the experimental data consists of measurements made by NA49 collaboration in central Pb-Pb collisions 
at beam momenta of 30, 40, 80 and 158A GeV, corresponding to -s/sjvat = '^■6' ^-^i 12.3 and 17.2 GeV respectively. 
The acceptance region in rapidity and transverse momentum covers a typical range from midrapidity to projectile 
rapidity and from to 1.5 GeV/c respectively. The overall hadron multiplicities, quoted in referenced papers, were 
obtained using forward-backward symmetry in rapidity and by extrapolating the yields to full phase space. All results 
were corrected for the feed-down from weak decays, e.g. tt" multiplicity does not include pions produced in decays of 
A hyperons and Kg mesons. 

Gentral collisions were selected by a trigger using information from a downstream calorimeter, which measured the 
energy of the projectile spectator nucleons. Whilst at 30, 40 and 80A GeV all pubhshed results refer to the 7.2% most 
central collision sample, at 158A GeV different centrality selections (5%, 10% and 20% most central collisions) were 
used to measure various hadronic species. In this analysis we have rescaled all published multiplicities at 158 A GeV 
to the corresponding ones at 5% most central collisions assuming that for the considered central collisions the hadron 
yield is proportional to the mean number of participant nucleona. The resulting scaling factors are 1.08 and 1.32 for 
10% and 20% most central collisions respectively [lll2^. 

As far as AGS data at 11. 6A GeV is concerned, we have used both multiplicities measured by the experiments 
and extrapolations of measured rapidity distributions made in ref. 9] at 3% top centrality. For A we have made a 
weighted average of the multiplicities measured both at 5% top centrality by E896 and E891 For the former, 
the quoted experimental error was only statistical so that we have added a 10% systematic error resulting in a value of 
16.7±0.5±1.7. For the latter, we have used the extrapolated value in ref. i9|| of 20.34±2.74. The error on the weighted 
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average has been rescaled by 1.25 (i.e. \/x^) according to the PDG weighting method in case of discrepancy between 
different measurements The obtained average has been rescaled by a factor 1.02 to convert it from 5% to 3% top 
centrahty by assuming a linear dependence on the number of participants and by using the tables in ref. [2^ . Since 
the A to A ratio has been measured only at midrapidity ^2d\ , we have obtained a A 47r multiplicity assuming that the 
double ratio ((A)/(A))j,=o/((A)/(A)) is the same at SPS and AGS energies. The final experimental multiplicities and 
ratios used in our analysis are shown in tables ^ lllll IIVI and ]V\ 

In order to test the effect of the cut in rapidity on the resulting statistical model parameters (discussed in detail 
in Sect. 4) we have also determined the yields integrated over limited {Ay = 1 and Ay = 2) rapidity windows around 
midrapidity in Pb-Pb collisions at 158A GeV. This has been done by fitting the rapidity distributions measured by 
NA49 to a Gaussian or the sum of two Gaussians, with area and width as a free parameters and central values set 
to zero. The results are shown in table IVII The quality of the fits is quite good, except for pions due to a couple of 
points near midrapidity; yet, this discrepancy does not affect significantly the integrated yield. In fact, it must be 
stressed that the main goal of these fits is to estimate an integral and not to reproduce accurately the shape of the 
distributions over the full measured range. We have also checked that the extrapolations to full phase space are in 
good agreement with published measurements. 



IV. ANALYSIS RESULTS 



The analysis has been carried out by looking for the minima of the ■ 

X^=E ^' ^ (4) 

i 

where is the multiplicity of the i"^ measured hadronic species and (7^ is the sum in quadrature of statistical and 
systematic experimental error. 

The theoretical multiplicities are calculated according to Eq. 10) with the decay chain stopped to match the 
experimental definition of multiplicity to properly compare theoretical and experimental values. This occurs in Pb-Pb 
collisions after electromagnetic and strong decays and before weak decays, whilst in Au-Au collisions at AGS the weak 
decays of A, S, S, O and Kg are included. 

The effect of the uncertainties on masses, widths and branching ratios of the involved hadrons on the fit parameters 
has been studied by the method described in ref. 0| and found to be negligible throughout. 

In order to cross-check our results and verify their robustness, we have performed the analysis with two independent 
numerical programs, henceforth referred to as A and B, which mainly differ with regard to the included resonances, 
their decay modes and branching ratios. 

The fitted parameters within the main scheme SIIM(75) are shown in table IVIII while the experimental and 
fitted multiplicities, along with the predicted yields of several hadron species are shown in tables HI ITU IIIII Hvl and 
IVl and figs, m m O and 21 Wc do not show any plot for the fit in Pb-Pb collisions at 30 A GeV because of the low 
number of data points. The quality of the fit is good throughout, as proved by the x^'s values and we do not see 
any clear discrepancy between data and model, with the remarkable exception of the A(1520) in Pb-Pb collisions at 
158A GeV. Due the 5a deviation from the statistical model prediction, the measured A(1520) yield has been removed 
from the fitted data sample as it could have biased the fit itself. We argue that this disagreement owes to its short 
lifetime (F — 15.6 MeV) compared with all other used particles. If the kinetic freeze-out occurs after some suitable 
delay from the chemical freeze-out, one can indeed justify the low measured A(1520) yield as the effect of the elastic 
reinteractions of its decay products. The quality of the fits is further illustrated in fig. |S| where the measured and 
fitted ratios (7r~)/(Np)and (K"'')/(7r+)are plotted as a function of ^snn', these ratios have been chosen as it has been 
goposed that their energy depedence plays an important role in the search for deconfinement onset at SPS energies 

The observed differences in the fit parameters between A and B are of the order of the fit errors. They may be 
considered as an estimate of the systematic error due to uncertainties in the implementation of the model. The first 
set of parameters in table IVIII have been obtained by using the full data, whilst the second set has been obtained 
by using the maximal common set of particles measured in the four collisions, that is tt"*", K"*", K~, A, A and the 
participant nucleons (net baryon number) B. By comparing fit results in the same analysis (A with A and B with 
B), it can be seen that the obtained parameters are in good agreement and only in one case a discrepancy larger than 
one standard deviation (7s in Pb-Pb at 158^ GeV) is observed; this demonstrates the robustness of the results. 

We have also included in tables m IIIII IIVI and the prediction for the yield of the recently discovered 0+ 
pentaquark baryon (uudds) by using as input mass m — 1540 MeV and J = 1/2. According to the SIIM(75) model. 
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in the Boltzmann limit this simply reads: 

(ne+) = ^m^TKa (^) eMt^s/T + iiq/T + ^ig/T] (5) 

if we disregard feeding from possible excited states. 

As the number of data points in Pb-Pb collisions at 30^ GeV was not sufficient to determine the four free parameters 
unambiguously, we have forced T to lie on the parabolic chemical freeze-out curve (|19|l in Sect. 5, interpolating the 
other four points in the — T plane. This method has proved to be able to provide unambiguous solutions for the 
remaining three free parameters. 

The first quoted error beside the best-fit value in table lVIll is the error coming out from the fitting program (inferred 
from the analysis of the = Xmin + 1 l^'^^l contours) whereas the second error is the fit error rescaled by a factor 
•\/x^~/do/ where dof is the number of degrees of freedom. We deem that the latter is a more realistic uncertainty on 
the parameters with respect to the fit error because of the "imperfect" x^^^/dof values, expected to be 1 on average 
if the model correctly matched physical reality. The argument, which is the same used in the Particle Data Book |22| 
when averaging discrepant data, is as follows: if Xmin/^o/ ^ 1, then the model cannot reproduce the data at the level 
of accuracy relevant to the experimental errors; on the other hand, this would be the case if experimental errors were 
larger and, particularly, if they were rescaled by a common factor S so that: 

i 

With this simple rescaling of the x^, the best-fit parameters would be unchanged, whereas their relevant errors would 
scale up by a factor S. In fact, the new covariance matrix C for the parameter vector X is related to the x^ around 
the minimum through: 

X''{X) = x^in + XofC-\x - Xo) (7) 

Also: 

X^{X) = xLn + iX~XofC-\X-Xo) (8) 

and, being x^ = X^/S'^i one finds C = S^C. Therefore, since this covariance matrix would be the outcome of a 
standard quality fit, with xtxin/'^^f ~ 1' the errors rescaled by S can be regarded as the sensible minimal uncertainty 
on the parameters. These rescaled errors have indeed been used in all of the plots in this paper. 

A major issue in the multiplicity fits is where to stop the inclusion of heavy light-flavoured resonances contributing 
to measured particle yields in Eq. Q . The relevance of this cut-off is owing to the peculiar shape of the hadron mass 
spectrum, which rises almost exponentially between 1 and 1.7 GeV and drops thereafter probably due to the missing 
experimental information (see fig. ^ . Should the number of states keep on increasing exponentially, the problem is 
set of the physical meaning of the obtained parameters, which could be heavily affected by the ignorance of further 
hadronic states. In fact, although the production of resonances decreases exponentially with the mass, the effect on 
secondary light particles through the decay chain could be balanced and even exceeded by the increasing number of 
states. We have thus checked the stability of the obtained parameters in the four collisions by varying the cut-off on 
the mass spectrum in a range where we are reasonably confident on the complete experimental knowledge and the 
number of states apparently rises exponentially (i.e. up to 1.7-1.8 GeV) and repeating the fit. As shown in fig. [7| 
the fitted T, fj,B and 75 in Pb-Pb at 158^ GeV are indeed fairly constant from 1.3 to 1.9 GeV. Furthermore, the 
outcoming primary yields of some measured particles tend to saturate at cut-off masses of about 1.8 GeV, implying 
that the contribution of resonance decays to secondaries (needed to keep the final multiplicity close to the measured 
value) settles down as well and the inclusion of heavier states yields a more and more negligible contribution. This is 
a clear indication of the significance of the fit results. A similar pattern occurs in all other examined collisions. 

A major result of these fits is that 75 is significantly smaller than 1 in almost all cases (with a possible exception 
at 30 A GeV, see table IVII|I . that is strangeness seems to be under-saturated with respect to a completely chemically 
equilibrated hadron gas. This confirms previous findings 0, |2^ . There is a considerable interest and ongoing 
investigations about this deviation of the data from the fully equilibrated hadron gas, particularly motivated by the 
fact that strangeness production is considered as a possible QGP signature. It is therefore worth to examine and test, 
with the presently available large data sample, different scenarios which have been put forward to account for the 
under-saturation of strangeness. 
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A. Full equilibrium and midrapidity ratios 



As we have seen, fits to full phase space multiplicities within the SHM yield 7s < f in most cases. However, 
good tests of the same model without extra strangeness suppression (i.e. assuming 75 — 1) have been obtained by 
fitting ratios of hadronic yields within a limited rapidity range around midrapidity at top SPS energy Q- This is an 
appropriate method of estimating the parameters of the sources only if the boost-invariant Bjorken scenario holds, 
at least as a good approximation, over a large rapidity interval (say Ay ~ 6) because, in this case, particle ratios at 
midrapidity are the same as in full phase space. However, rapidity distributions of hadrons at SPS energies do not 
feature boost-invariance P, Is^l and a cut at midrapidity can artificially enhance heavy particle yield 

with respect to light ones (see fig. (Sj), as long as their kinetic frcczc-out occurs at the same temperature and the 
leading baryon effect can be neglected. In the statistical model of a single fireball this can be easily understood, for 
the width of the rapidity distribution decreases as a function of mass according to (in the Boltzmann approximation): 



diV / . 2mT^ 2T3 \ ^ , , , 

— cx m-^T H \ ^ exp -TOCOshy/T (9) 

dy \ coshy cosh"^ y / 

Yet, it is worth testing the effect of the rapidity cut on measured distributions rather than using arguments based 
solely on the statistical model. Therefore, we have fitted, within the scheme A, the integrated yields over limited 
rapidity windows measured in central Pb+Pb collisions at 158A GeV by NA49 (see Sect. 3 for details) and quoted 
in table I VI III as well as the yields measured by WA97 [sj over a Ay = 1 window around midrapidity. We first note 
that, according to table I VI Hi the integrated yields over Ay = 1 measured by NA49 and WA97 are in good agreement 
with each other. Since the fit to the statistical model gave 75 ~ 1 for WA97 data, the same is expected for the 
integrated NA49 yields over the same rapidity window. This is indeed what we find, as shown in table Hxl While 
temperature and baryon-chemical potential are essentially unchanged, the best-fit value of 75 is closer to 1 than that 
obtained in full phase space in fit A (see table l|VII|) and it is also compatible with 1 within the error. 

We then conclude that the superfiuity of 7^ in analysis of midrapidity particle yields, at least at top SPS energy, is 
likely to owe to the artificial enhancement of strange particles with respect to lighter non-strange ones, induced by the 
cut on rapidity. The fact that 75 ~ 1 for midrapidity yields is then not an indication of a fully equilibrated hadron gas 
at midrapidity; even if such equilibrated fireball existed at the estimated kinetic freeze-out temperature of T « 125 
MeV 33], the Ay — 1 window would be too narrow for a correct estimation of chemical freeze-out parameters (see 
fig. IHl because lighter particles would be cut down significantly. 



B. Strangeness correlation volume 



To account for the observed under-saturation of strangeness, a picture has been put forward in which strangeness is 
supposed to be exactly vanishing over distances less than those implied by the overall volume V 34] . We henceforth 
refer to this version of the statistical model as SHM(SCV). Following the description of the model in Sect. 2, this 
means that the produced clusters or fireballs emerge with S = and they are not allowed to share non-vanishing 
net strangeness. Assuming, for sake of simplicity, that all clusters have the same typical volume Vc and that the 
equivalence of the set of clusters to a global fireball still applies for baryon number and electric charge (but not to 
strangeness) the following expression of the primary average multiplicities can be obtained (see Appendix B): 



/ iTlr+'eM-n^f^^,/T + n^^nB,/T + n^^QQ,/T]^-^ (10) 



where 



Zr = ^ f dcj, cxp &±il^ j d^p log(l ± e-V^^/^+'^-^./^+'^«Q./^-'^^^.)±i] (11) 



27r 



is the so-called strange canonical partition function of a single cluster. As usual, in the above equations, the upper 
sign is for fermions and the lower for bosons. 

If Vc is sufficiently small, the multiplicities of strange hadrons turn out to be significantly suppressed with respect 
to the corresponding grand-canonical ones due to the enforcement of exact strangeness conservation in a finite system, 
an effect called canonical suppression. Furthermore, the suppression features hierarchy in strangeness, namely it is 
stronger for (S* = 3) and S (S* = 2) than for kaons and A's, so it can be argued that this can account for the actually 
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observed hierarchical pattern of extra strangeness suppression which goes hke 7^ for open strange particles. The 
discriminating difference between this picture and our main scheme SHM(75) described in Sect. 2 is concerned with 
hidden strange particles such as 4>, which do not suffer canonical suppression, so that its theoretical multiplicity in 
SHM(SCV) turns out to be simply the same as in a cluster with volume V, that is given by the formula ^ without 
7I suppression. 

We have made a test of this model by fitting the data sample of full phase space multiplicities in Pb-Pb collisions at 
158 A GeV fixing 75 = 1 and determining the parameters T, V, hb and / = Vc/V within the scheme A. The results 
are shown in tables Ixl and IVl The quality of the fit is worse with respect to the SHM(75') model mainly because of 
the underestimated pion yield and the larger of 4>. The latter is expected, as has been mentioned. As far as pion 
discrepancy is concerned, the deviation stems from the very fact that they are the only non-strange particles in the 
fit. The minimization procedure tries to accommodate the relative ratios among strange hadrons by fixing Vc and T, 
then it tries to set the overall normalization V and at this stage a competition sets in between the set of strange and 
non-strange particles. Since pions are only two entries, the fit prefers to keep them low rather than raising the whole 
set of strange particles. 

Our result suggests that, for the local strangeness correlation to be an effective mechanism, the cluster volume 
should be of the order of 2.5% of the overall volume. Otherwise stated, strange quarks should have not propagated 
beyond a distance of about 4 fm from the production point up to chemical freeze-out, if we take the overall volume 
of about 3 10'^ fm"^ as coming out from this fit where hadrons are pointlike particles. 



C. Superposition of NN collisions with a fully equilibrated fireball 

In this picture, henceforth referred to as SHM(TC), the observed hadron production is approximately the superpo- 
sition of two components (TC): one originated from one large fireball at complete chemical equilibrium at freeze-out, 
with 75 — 1, and another component from single nucleon-nucleon collisions. In fact, according to simulations based 
on transport models, a significant fraction of beam nucleons interacts only once with target nucleons |35l |. With 
the simplifying assumption of disregarding subsequent inelastic collisions of particles produced in those primary NN 
collisions, the overall hadron multiplicity can be written then as: 

{nj)^{N,)(n,)NN + {nj)v (12) 

where {rij) jvat is the average multiplicity of the j*'' hadron in a single NN collision, (A^c) is the mean number of single 
NN collisions giving rise to non-re-interacting particles and {nj)v is the average multiplicity of hadrons emitted from 
the equilibrated fireball, as in Eq. with 75 — 1. The {nj)^^ term can be written in turn as: 

Z2 ^ 2Z{A-Z), . 
{nj)NN = -^{nj)pp H ^2 m)nn H ^2 m)np (13) 

Since it is known that in NN collisions strangeness is strongly suppressed ths idea is to ascribe the observed 
under-saturation of strangeness in heavy ion collisions to the NN component, leaving the central fireball at complete 
equilibrium, i.e. with 75 = 1. Of course, this is possible provided that (Ac) is sufficiently large. This production 
mechanism has probably some consequences on the final rapidity and momentum distributions of the different species, 
whose calculation goes certainly beyond the scope of this paper. Instead, we have confined ourselves to integrated 
multiplicities and tried to fit T, V, jiB of the central fireball and (Ac) by using NA49 data in Pb-Pb collisions at 
158 A GeV within the scheme A. 

To calculate {nj)NN we have used the statistical model and fitted pp full phase space multiplicities measured at 
^/s = 17.2 GeV (i.e. the same beam energy) by the same NA49 experiment. For np and nn collisions, the parameters 
of the statistical model determined in pp are retained and the initial quantum numbers are changed accordingly. 
Theoretical multiplicities have been calculated in the canonical ensemble, which is described in detail in ref. [l4l| . 
Instead of the usual 75 parametrization, the new parametrization described in ref. |Q| has been used in which one 
assumes that some number of Bpairs, poissonianly distributed, hadronizes; the extra strangeness suppression parameter 
75' is thus replaced by the mean number of these Bpairs, (ss). 

The results of the fit are shown in table lXlj along with fitted and predicted hadron multiplicities, including the 8(1540) 
pentaquark baryon, and in fig. El The temperature value is significantly higher than in pp and pp collisions at higher 
energy, an effect already observed for center-of-mass energies below 20 GeV dM- We conjecture that this is a 
possible indication of a beginning inadequacy of the canonical ensemble at low energy, where exact conservation of 
energy and momentum should start to play a significant role. Perhaps this is the point where the microcanonical 
hadronization of each cluster is a more appropriate approach. 

The results of the fit to Pb-Pb collisions are shown in table |3 The fit quality, as well as the obtained values of T, 
/is, are comparable to the main fit within the SIIM(75) model. The predicted number of "single" NN collisions is 
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about 50 with a 16% uncertainty. Thus, only 260 nucleons out of 360 contribute to the formation of large equilibrated 
fireballs. The percentage of primary hadrons stemming from NN collisions varies from 14% for pions to 27% of p's and 
protons and to 0.5% of 17's. It should be pointed out that the fitted parameters are affected by a further systematic 
error owing to the uncertainty on the parameters of the statistical model in NN collisions, which are used as an input 
in the Pb-Pb fit. However, because of exceeding computing time needs, it has not been possible to assess these errors. 

In a simple geometrical picture, the single-interacting nucleons are located in the outer corona of the portion of 
colliding nuclei corresponding to the observed number of participants. As the projected (on the collision's transverse 
plane) radial nucleon density is: 

^ = 4:7rr\/ R'^ - r^no (14) 
dr 

where uq = 0.16 fm^'^ is the nucleon density and R ~ 6.45 fm is the radius of the portion of colliding nucleus 
corresponding to a participant number of 180, the 50 single-interacting nucleons should lie between 4.84 and 6.45 
fm. This simple estimate is in approximate agreement with the calculations performed with the Glauber model |35| | . 
Since the number of single-interacting nucleons is expected to be weakly dependent on center-of-mass energy, the fits 
to this two-component model should yield consistent values of (Nc) at 30, 40 and 80^4 GeV collisions. However, no 
measurement of hadron production in NN collisions at those energies is available and this question cannot be tackled 
for the present. 



D. Non-equilibrium of hadrons with light quarks 



An extension of the statistical model has been proposed where QGP hadronization is essentially a statistical 
coalescence of quarks occurring at an energy density value which does not correspond to a hadron gas at equilibrium 
|36j |. In this model two non-equilibrium parameters are introduced for the different types of quarks, for u, d quarks 
and 7s for strange quarks (the difference between 7^ and 7s is explained below). The multiplicity of each hadron thus 
reads: 



(2J,- + 1)^ 
(27r)3 



00 „ 

E 7r=7r' J d'p exp[-n^p2 + m2/T + n/x • q,/T] (15) 



where Us is the number of valence s quarks and Uq the number of valence u, d quarks; fi and qj are as in Eq. By 
defining: 



75 = - (16) 



the Boltzmann limit of average multiplicity reads 

, , _ (2 J, + l)V 



(27r)3 



-75=7^^1 / d^P exphy'p^+m^/T + /x • q,/r] (17) 



where Bj is the baryon number, as long as mesons have two and baryons have three valence quarks. By comparing this 
formula with the Boltzmann limit of Eq. |^ it can be realized that the introduction of a light-quark non-equilibrium 
parameter amounts to introduce in the statistical model an overall enhancement (or suppression) of baryons with 
respect to mesons, unlike in the model SHM(75). We henceforth refer to this model as SHM(757q). 

The parameter 7^ has a definite physical bound for bosons which can be obtained by requiring the convergence of 

the series 'Y^'^=i^ilq''^)^^v{'^^^lT + Nfj, ■ (ij/T) for any value of the energy. If the number of u, d quarks to be 
hadronized is so large that 7g is to attain its bounding value, a Bose condensation of particles in the lowest momentum 
state sets in. For low strangeness and electrical chemical potentials, such as those found in the present analysis, the 
bounding value is 7, — exp(m/2T) where m is the neutral pion mass, e.g. jq ~ 1.5 for T ~ 160 MeV. 

With the introduction of jq as an additional free parameter, there are 5 parameters to be determined in the model. 
This makes the minimization procedure rather unstable because it becomes easier to be trapped in local minima. To 
avoid this, we have performed 4 parameter fits with fixed values of 7^ varying from 0.6 to 1.7 in steps of 0.1. This 
method allows a clearcut determination of the absolute minimum. 

The results of these fits are shown in fig. ^]in terms of the minimum obtained for fixed jq. The round dots 
show the minimum x^'s achieved by using the main sample of multiplicities in Pb-Pb collisions at 158 A GeV. We find 
a steady decreasing trend in the value of best-fit temperatures varying from ~ 187 MeV at 7, = 0.6 to ~ 140 MeV at 
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7q = 1-6. The fitted temperature at 7^ = 0.5 reaches the upper Hmit of 200 MeV, which is the maximum allowed in 
the model to prevent from being critically dependent on the hadron mass spectrum cut-off. On the other hand, the 
best- fit values of hb/T and 75 are rather stable and about the same found in the main fit with 7^ = 1. The number 
of terms in the series (I15f) has been truncated to 5 for all particles; the contribution of further terms has been found 
to be neghgible throughout. 

It is seen that the absolute minimum falls in the region of pion condensation, marked by a vertical line at 
7g ~ 1.62, with — 13 and T ~ 140 MeV. This finding is in agreement with what is found in ref. j36i |. However, 
there is also a local minimum at the lower edge 7^ = 0.6, with a temperature of 187 MeV, which is only one unit of 
higher than the absolute minimum. This indicates that the absolute minimum could be rather unstable against 
variations of the input data and this is in fact what we find by varying down the pion multiplicities by only 1 a. For 
this case, the minimum x^s are shown in fig. 1101 as triangular dots and the absolute minimum now lies at 7^ = 0.6 
instead of at the pion condensation point. 

In view of the instability of the fit, and of the small relative improvement in comparison with the main fit, 
we conclude that there is so far no evidence for the need of this further non-equilibrium parameter. However, it is 
interesting to note that this model predicts an enhanced (if 7^ > 1) or suppressed (if 7^ < 1) production of the 0+ 
pentaquark baryon with respect to the other versions of the statistical model, even though they agree in reproducing 
the multiplicities of all other hadrons. This is owing to an additional 7^ factor for this special hadron having five 
valence quarks. From Eqs. (|15I16|I one gets, in the Boltzmann limit: 

'Y"^a/c'\/ /JJl\ 

(ne+) = ^^^m^TKa (- j exp[fiB/T + mq/T + fis/T] (18) 

In fact, as can be seen in table^the predicted yield of 8^" at the global minimum 7, ~ 1.62 is more than a factor 2 
higher than for SHM(75). 

As a final remark, we stress that minimum x^ fits are very useful tools to get information on the state of the source 
at chemical freeze-out, but, as already emphasized in Sect. 3, the simple multiplicity analysis with global parameters 
resides on an idealization of the collision (e.g. the assumed existence of an EGG) which cannot exactly fit physical 
reality and discrepancies are to be expected anyway. Thus, a new mechanism or of a modification of the basic scheme 
proves to be relevant only if it leads to a major improvement of the agreement with the data. Slight improvements of 
the x^j whenever their significance is beyond its expected statistical fluctuations, cannot be seriously taken as a proof 
of the validity of a particular scheme. 

V. ENERGY DEPENDENCE 

The statistical model does not make any prediction on the energy dependence of hadron production; its relevant 
parameters have to be determined separately for each energy and reaction type. Nevertheless, the analysis of the data 
within this model may help in the study of energy dependence of hadron production because it effectively reduces the 
full experimental information on numerous hadron yields to only few parameters. Glearly this reduction should be 
taken with care, were not for the approximate validity of some relevant assumptions, such as the reducibility to EGG 
(see Sect. 2). Furthermore, the reduction procedure may remove or dilute essential physical information. With these 
caveats in mind, in this section we discuss the energy dependence of the chemical freeze-out parameters extracted 
from the data (by using the full data sets) within our main SHM(75) approach. 

The chemical freeze-out points in the /is — T plane are shown in fig. The RHIG point at VsjvAf ~ ^^O GeV, 
obtained fitting particle yield ratios at midrapidity, has been taken from ref. [ll|. The four points at beam energies 
of 11.6, 40, 80 and 158^ GeV have been fitted with a parabola: 

T = 0.167- 0.153Ai|, (19) 

where T and /is are in GeV. The Pb-Pb point at 30 A GeV has been forced to lie on the above curve, as has been 
mentioned in the previous section. The RHIG point calculated in ref. [llj | is in good agreement with the extrapolation 
of the curve H19() . 

In the search for deconfinement phase transition, strangeness production is generally believed to be a major item, 
especially if an anomalous abrupt change was found as a function of centre-of-mass energy or other related quantities. 
A possible indication of it in Pb-Pb collisions at the low SPS energies was reported on the basis of the observed 
energy dependence of several observables 0|. Particularly, the (K"'")/(7r+) ratio shows a peaked maximum at about 
30A GeV. One may expect that this anomaly should be reflected in the energy dependence of 75 parameter fitted 
within SHM(75) scheme. This dependence is plotted in fig.ll2land in fact a maximum shows up at 30 A GeV. Although 
the error bars are large enough so as to make 75 seemingly consistent with a constant as a function of centre-of-mass 
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energy, it is important to note that the dominant systematic errors on experimental data at SPS energies are essentially 
common. Therefore, the errors on the model parameters at different SPS energies turn out to be strongly correlated, 
hence fitted 75 's are expected to move up or down together. 

In order to further study strangeness production features, we have also compared the the measured (K+) / (tt^) ratio 
(including the preliminary RHIC result at Vs^^ — 200 GeV ^3) with the theoretical values in a hadron gas along 
the freeze-out curve (|19f) as a function of the fitted baryon-chemical potential for different values of 75 (see fig. I13II . 
The calculated dependence of (K+)/(7r+)on p^B is non-monotonic with a broad maximum at [Lb — 400 MeV (i.e. 
Ebeam — 30 A GcV) 38] . Taking into account that systematic errors at different energies in Pb-Pb collisions are fully 
correlated, we can conclude that the data points seem not to follow the constant 75 lines. 

In fact, the anomalous increase of relative strangeness production at 30 A GeV can be seen also in the Wroblewski 
variable A5 — 2(ss)/((uii) + (dd)), the estimated ratio of newly produced strange quarks to u, d quarks at primary 
hadron level, shown in fig. 1141 and table IVIll The calculation of newly produced quark pairs is performed by using 
the statistical model best fit values of the various hadron multiplicities, so the obtained A5 values are somehow 
model-dependent. Nevertheless, this variable features a very similar behaviour as the ratio (K+)/(7r+)and attains a 
maximum value of 0.61 at 30A GeV, very close to that predicted for 75 = 1. 

These deviations from a smooth behaviour of strangeness production are certainly intriguing, yet the analysis within 
the SHM will be more conclusive in this regard with a larger data sample at 30A GeV and at the forthcoming 20^1 
GeV data. 



VI. SUMMARY AND CONCLUSIONS 

We presented a detailed study of chemical freeze-out conditions in ultra-relativistic heavy ion collisions at projectile 
momenta of 11.6 (Au-Au at AGS), 30, 40, 80 and 158A (Pb-Pb at SPS) GeV, corresponding to nucleon-nucleon centre 
of mass energies of 4.8, 7.6, 8.8, 12.3, 17.2 GeV respectively. By analyzing hadronic multiplicities measured in full 
phase space within the statistical hadronization model, we have tested and compared different versions of this model, 
with special emphasis on possible explanations of the observed strangeness under-saturation at the hadron level. 

It is found that version of the model referred to as SHM (75), where a non-equilibrium population of hadron carrying 
strange valence quarks is allowed, fits all the data analyzed in this paper. We have also shown that the seeming full 
chemical equilibrium found for central Pb-|-Pb collisions at 158^1 GeV by using particle yields integrated over a limited 
region around mid-rapidity is most likely an artefact of the kinematical cut. 

We have tested a model (SHM(TC)) in which hadron production is pictured as stemming from two independent 
components: a fireball (or a set of fireballs) at full chemical equilibrium and single nucleon-nucleon interactions. This 
model can fit the data at 158A GeV if the number of collisions is around 50 with a sizeable uncertainty. So far, it 
cannot be confirmed at other energies due to the lack of the precise data on NN interactions. 

A model in which strangeness is assumed to vanish locally [SJI yields a worse fit to the data with respect to SHM(75) 
and SHM(TC). 

Finally, we have also tested a model in which it is allowed a non-equilibrium population of hadrons carrying both 
strange and light valence quarks. We have found the present set of available data does not allow to establish whether 
a further non-equilibrium parameter is indeed needed to account for the observed hadron production pattern. A 
discriminating prediction of this model with respect to the SHM(75') is an enhanced production of the recently 
discovered pentaquark baryon due to the additional factor 7^. 

Energy dependence of chemical freeze-out parameters has been discussed based on the results obtained with our 
main version of the model SHM (75). The evolution of the freeze-out temperature and baryon-chemical potential 
is found to be smooth in the AGS-SPS-RHIC energy range. The strangeness-suppression parameter 7c is found to 
smaller than one (7s ~ 0.8) for most of the studied collisions which confirms previous findings |E0ill3i with an 
indication of a maximum at 30A GeV, where 7s is found to be close to one. The significance of this maximum is related 
to the correlation between errors on hadron yield measurements at different energies. The interpolated dependence 
of relative strangeness production on energy and baryon-chemical potential, as measured by the (K"'")/(7r+)ratio 
and the Wroblewski factor Ag, features a broad maximum at about 30 A of beam energy. However, the experimental 
measurement of (K+) /(7r+)and the estimated As value in Pb-Pb collisions at this energy seemingly exceed the expected 
values for a fixed js- 
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APPENDIX A - FROM CANONICAL TO GRAND-CANONICAL ENSEMBLE 

The canonical partition function of the z'^ cluster can be written as a multiple integral over the interval [— tt, tt] 

^i(QO = I d3<^e'Q' <^exp[^^((/,)] (20) 

where Qi — {Qi, Bi, Si) is a vector having as components the electric charge, baryon number and strangeness of the 
cluster, (f> = {(j)Q,(l>B,4's) and F{<p) reads: 



(21) 



Vi is the volume and T the temperature of the cluster; the sum runs over all hadronic species j and the charge 
vector of the j^^ hadron; the upper sign applies to fermions, the lower to bosons. The probability distribution required 
for the reduction to EGC to apply reads |a, ^3 ■ 

r> ^ n»^»(Q»)'^i:.Q.,Q 

Z^Qi,....Q„ Hi Q.,Q 



In this case, the overall multiplicity of the j'^ hadron is given by [T3l| : 

dX 



= AiogZ(Q)^^^ (23) 



where: 



where Z{Q,) is the canonical partition function of the equivalent global cluster: 

Qi,...,Qjv i 

and Xj is a fictitious fugacity parameter. Formally, this turns out to be the same function as in Eq. (|20|) with V = Vi 
replacing Vi and Q = replacing Q^. If ^ is large, the canonical partition function can be approximated by the 

leading term of an asymptotic saddle point expansion. Setting exp(— i0/j) = Zk, the canonical partition function can 
be written as: 

^(Q) = (2^ [ n / ^^1 n (25) 

^ k=i ^ k=\ 

= E ^'"'llT [ 1"S(1 ± e-^/^/^ z^)^^ (26) 
The saddle-point expansion is carried out by requiring the logarithmic derivative of the integrand to vanish: 

-^ + ^=0 fc = 1,2,3 (27) 
Zk ozk 

The solutions of this equation are indeed the grand-canonical fugacities Xk = exp(/Zfe/T). The function Fc{\) coincides 
with the logarithm of the grand-canonical partition function log Zgc, therefore the equation (|77l) expresses the equality 
between the average charge in the grand-canonical ensemble Xkdlog Zgc/dXk and the initial value Qk- The canonical 
partition function now becomes, at the second order of the expansion: 

3 3 

Z{Q) ^ exp[F,(A)] [ n A,-^'=] [fl <f —] cxph(z - A)^H (z - A)/2] (28) 

k=i ^ k=i-' 

where H is the Hessian matrix in Zk = Xk- The first exponential factor is just the grand-canonical partition function 
Zgc calculated for the fugacities . The average multiplicity of the j*** hadron species can now be calculated by using 
Eq. (|23|) taking the approximated expression (|28|) of the canonical partition function. Retaining only the dominant 
contribution, one just obtains the grand-canonical expression of the average multiplicity as expressed in Eq. (Q. 
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APPENDIX B - PROOF OF EQUATION JTHl) 



The argument closely follows the previous one. The main difference is the request of vanishing strangeness for each 
cluster. Thus, the configurational probabilities (|22|l w turn to: 



and the average multiplicity of the j*"^ hadron species now reads: 



d 

' Ql,---,QiV i 



A,=l 



where Xj is a fictitious fugacity. Let us now work out the expression: 



N 



Qi,...,Qn 1=1 



(29) 



(30) 



(31) 



by assuming that all the clusters have the same volume Vc and temperature T. This can be done rewriting the canonical 
partition functions Zj like in Eq. H2U|) . using the integral representation of the Kronecker's delta and eliminating the 
redundant strangeness conservation constraints. Thus, expanding the vector Q in its components, the Eq. H31|) 
becomes: 



E 

B\,Q\,...,B]\j ,Q 
N 

xn 



i=l 



27r 



d0is r d(/)jQ 



27r 



27r 



27r 



27r 



exp[f^((/.,B,0»Q,0zs)] 



(32) 



where F is the function in Eq. (|21|l with T4 replacing Vi. Note that this function is the same for all clusters, being T 
and Vc constant. We can now carry out the sum over all the integers Qi in Eq. (|32|l and get: 



E 



N 

e-iE.s.(,/,B-0.i,)-iE.Q.WQ-0.Q) ^ J|(27r)3,5(0B 

Bi,Qi,...,Bn:Qn i=l 

so that the integration over and in Eq. H32|l can be easily done and one is left with 

N 



(33) 



C 



As the function F is the same for all clusters, this can be written also as: 



c = 



2ti 



2ti 



e.S0.+.Q0«^ exp[F(</)s,. 



V/V, 



(34) 



(35) 



being V = Vi — NVc- For large volumes, we can approximate C by means of the saddle-point expansion of the 
integrals over (pB and 0q like in Appendix A. Thus, similarly to Eq. (^5)1 : 



I j\ If [E ^^7^ / ± Af^ Ag^e-V^/--^^-0^^ 



The function: 



^ J[ ^ exp [E ^^^1^ / d3p log(l ± Af^ Ag^e-v^/--'^^^0^^ 



(36) 



(37) 



is defined as the strange canonical partition function. The multiplicity of the hadron j can now be calculated by 
means of Eq. Ip}|) by using Eqs. and the definition 133). What is obtained is just Eq. ((117)) 
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TABLE I: Comparison between measured and fitted particle multiplicities, in the framework of SHM(75) model in central 
Au-Au collisions (3%) at a beam energy of 11. 6^1 GeV. Also shown the predicted multiplicities of the main hadron species. 





Reference 


Measurement 


Fit A 


Fit B 


Np 


[39] 


363 ± 10 


361.7 


360.6 


1 _|_ 




1.23 ± 0.13 


1.277 


1.224 


4- 


40] 


133.7 ± 9.93 


134.9 


140.0 








176.9 


182.2 


n 






163.5 


163.2 


K+ 


[39] 


23.7 ± 2.86 


18.80 


18.81 


K" 


[39] 


3.76 ± 0.47 


3.890 


3.539 








11.68 


11.68 








8.073 


6.340 


LO 






4.870 


3.659 








0.3287 


0.3489 


v 






0.2997 


0.2437 


_l_ 






7.707 


10.39 


p~ 






9.164 


12.55 








8.517 


11.59 


K*+ 






3.555 


3.512 


K*~ 






0.6179 


0.5145 


K " 






3.766 


3.801 


K*° 






0.5555 


0.4628 


P 






172.2 


171.4 


P 






0.02851 


0.02465 


A++ 






25.39 


24.51 


A 






0.003071 


0.00222 


A 


[23, J4] see text 


18.1 ± 1.9 


19.82 


20.71 


A 




0.017 ±0.005 


0.01601 


0.01645 


E+ 






4.840 


4.784 


E" 






5.457 


5.453 


EO 






5.163 


5.106 


E" 






0.003445 


0.00321 


E+ 






0.002793 


0.00259 


E° 






0.003115 


0.00288 








0.56067 


0.5564 








0.54670 


0.5387 








0.002133 


0.00248 








0.002392 


0.00280 


n 






0.01352 


0.01459 


n 






0.0003569 


0.00056 


A(1520) 






0.7720 


0.6601 


e+(1540) 
G-(1540) 






1.86 
2.710~^ 


2.20 
1.8710"^ 
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TABLE II: Comparison between measured and fitted particle multiplicities, in the framework of SHM(7s) model, in central 
Pb-Pb collisions (7.2%) at a beam energy of 30A GeV. Also shown the predicted multiplicities of the main hadron species. 





Reference 


Measurement 


Fit A 


Fit B 


Np 


[2] 


349 ± 1 ± 5 


350.5 


350.5 






239 ± 0.7 ± 17 


228.4 


228.5 


TV~ 


[2] 


275 ±0.7 ±20 


256.5 


256.8 


n 






265.8 


251.9 


K+ 


m 


55.3 ± 1.6 ± 2.8 


49.83 


48.83 






16.1 ± 0.2 ± 0.8 


17.11 


20.72 








33.57 


34.85 


V 






23.74 


21.60 








15.45 


12.99 








2.571 


2.848 


v 






1.411 


1.341 


_l_ 






20.24 


22.68 


p~ 






23.09 


26.05 


A 

P° 






22.14 


25.11 


K + 






13.65 


13.45 


K*" 






4.006 


3.668 


K*" 






14.21 


14.20 


K " 






3.710 


3.386 


P 






138.0 


137.0 


P 






0.3650 


0.3803 


A -1-4- 






26.90 


25.23 


A 






0.07781 


0.07438 


A 






38.02 


40.25 


A 






0.3393 


0.3901 


I— .4- 






9.995 


9.935 


E" 






10.91 


10.89 


E" 






10.48 


10.40 


E" 






0.1005 


0.1106 


E+ 






0.08620 


0.09491 


E" 






0.09334 


0.1023 








2.422 


2.422 








2.378 


2.369 


=■ + 






0.07920 


0.09332 








0.08446 


0.09967 


n 






0.1587 


0.1799 








0.02067 


0.02959 


A(1520) 






2.167 


1.751 


e+(1540) 
9^(1540) 






2.84 
0.0018 


3.02 
0.0019 
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TABLE III: Comparison between measured and fitted particle multiplicities, in the framework of SHM(7s) model, in central 
Pb-Pb collisions (7.2%) at a beam energy of 40A GeV. Also shown the predicted multiplicities of the main hadron species. 





Reference 


Measurement 


Fit A 


Fit B 


Np 




349 ± 1 ± 5 


352.1 


351.6 






293 ± 3 ± 15 


285.5 


288.3 


TV~ 


fi] 


322 ± 3 ± 16 


314.7 


317.9 


n 

tt" 






330.4 


314.9 


K+ 


LJJ 


59.1 ± 1.9 ± 3 


51.22 


50.61 




LJJ 


19.2 ± 0.5 ± 1.0 


20.52 


20.33 


-rrC\ 






35.79 


36.24 


V 






30.48 


26.59 








23.00 


19.57 


t 


[22] 


2.57 ±0.10 


2.641 


2.644 


v 






1.858 


1.621 


_l_ 






28.84 


31.11 


p~ 






32.36 


35.01 


A 

P° 






31.27 


33.97 


K + 






15.55 


14.34 


K*" 






5.393 


4.860 


K*" 






16.06 


14.98 


K " 






5.045 


4.535 


P 






141.7 


141.9 


P 






0.9824 


0.9784 


A -1-4- 






29.05 


27.03 


A 






0.2181 


0.1969 


A 


[28] 


45.6 ± 3.4 


36.60 


37.36 


A 


[28] 


0.74 ± 0.06 


0.7223 


0.7297 


I— .4- 






9.655 


9.221 


E" 






10.40 


9.938 


E" 






10.05 


9.567 


E" 






0.2116 


0.2033 


E+ 






0.1853 


0.1783 


E" 






0.1985 


0.1901 








2.118 


1.948 








2.089 


1.917 


=■ + 






0.1285 


0.1207 








0.1358 


0.1277 


n 






0.1364 


0.1344 








0.02719 


0.02788 


A(1520) 






2.273 


1.688 


e+(1540) 
9^(1540) 






2.61 
0.0052 


2.32 
0.0045 
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TABLE IV: Comparison between measured and fitted particle multiplicities, in the framework of SHM(7s) model, in central 
Pb-Pb collisions (7.2%) at a beam energy of 80A GeV. Also shown the predicted multiplicities of the main hadron species. 





Reference 


Measurement 


Fit A 


Fit B 


Np 


fl] 


349 ± 1 ± 5 


352.0 


351.5 




m 


446 ± 5 ± 22 


420.3 


422.7 


TV~ 




474 ± 5 ± 23 


450.9 


453.7 


n 

tt" 






485.2 


457.6 


K+ 




76.9 ± 2 ± 4 


70.72 


69.87 




LJJ 


32.4 ± 0.6 ± 1.6 


35.96 


35.96 


-rrC\ 






52.80 


53.77 


V 






49.37 


43.47 








39.51 


34.43 


t 


[22] 


4.37 ± 0.14 


4.354 


4.353 


v 






3.196 


2.786 


_l_ 






47.42 


48.92 


p~ 






51.74 


53.54 


A 

P° 






50.79 


52.83 


K + 






22.81 


21.10 


K*" 






10.43 


9.600 


K*" 






23.28 


21.73 


K " 






9.892 


9.089 


P 






141.5 


142.5 


P 






3.379 


3.649 


A -1-4- 






30.07 


28.05 


A 






0.7623 


0.7439 


A 


[28] 


47.4 ± 3.7 


42.12 


42.85 


A 


[28] 


2.26 ± 0.35 


2.171 


2.328 


I— .4- 






11.23 


10.67 


E" 






11.82 


11.23 


E" 






11.56 


10.93 


E" 






0.6265 


0.6348 


E+ 






0.5643 


0.5729 


E" 






0.5961 


0.6024 








2.774 


2.505 








2.758 


2.485 


=■ + 






0.3279 


0.3154 








0.3428 


0.3299 


n 






0.2154 


0.2090 








0.06132 


0.06332 


A(1520) 






2.769 


2.028 


e+(1540) 
9^(1540) 






2.35 
0.022 


2.04 
0.021 
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TABLE V: Comparison between measured and fitted particle multiplicities, in the framework of various versions of the SHM, 
in central Pb-Pb collisions (5%) at a beam energy of 158^4 GeV. Also shown the predicted multiplicities of the main hadron 
species. The A(1520) multiplicity was not used in the fits (see text). For the SHM(7s,7q) model, the multiplicities are those 
calculated in the 7, — 1.6 fit. 



Particle 


Reference 


Measurement 


SHM(7s) fit A 


SHM(7s) fit B 


SHM(SCV) 


SHM(TC) 


SHM(7s,7,) 


Np 


HI 


362 ± 1 ± 5 


363.6 


363.7 


362.0 


364.2 


362.6 


7r+ 




619 ± 17 ±31 


551.5 


533.2 


502.7 


563.4 


578.7 


TV~ 




639 ± 17 ±31 


583.5 


565.6 


534.1 


595.3 


612.5 








638.2 


576.2 


585.4 


661.9 


661.6 


K+ 


fil 


103 ± 5 ± 5 


103.5 


103.9 


106.7 


99.98 


102.3 




fll 


51.9 ± 1.9 ±3 


59.57 


59.35 


59.54 


59.23 


57.77 




[28] 


81 ±4 


80.31 


81.13 


81.65 


78.19 


78.61 


77 






70.69 


62.90 


67.73 


76.94 


62.22 








54.93 


45.34 


48.75 


58.75 


43.01 


<l> 


[29J 


7.6 ± 1.1 


8.136 


8.676 


10.07 


9.088 


7.084 


V 






4.940 


4.461 


5.021 


5.471 


3.923 


P+ 






64.56 


61.67 


57.26 


67.18 


54.12 


p~ 






69.43 


66.82 


61.93 


72.29 


57.04 


p° 






68.86 


66.84 


61.44 


73.14 


57.00 


K*+ 






34.43 


32.05 


35.51 


32.63 


27.90 


K*" 






18.10 


16.47 


17.79 


17.61 


14.79 


K*o 






34.93 


32.81 


36.12 


33.07 


27.99 


K*o 






17.29 


15.69 


16.96 


16.84 


14.17 


p 






143.71 


142.9 


138.51 


147.7 


144.4 


p 






7.053 


6.877 


5.756 


7.721 


7.046 


A++ 






31.01 


28.32 


29.70 


30.70 


29.27 


A"- 






1.595 


1.393 


1.295 


1.716 


1.472 


A 


[28] 


53.0 ±5.0 


53.88 


56.22 


57.06 


49.53 


53.38 


A 


[28j 


4.64 ± 0.32 


4.976 


5.077 


4.698 


4.899 


4.878 


S+ 






14.45 


14.11 


15.28 


13.21 


14.57 


E" 






15.04 


14.71 


15.98 


13.73 


14.80 








14.78 


14.39 


15.67 


13.57 


14.73 


E" 






1.424 


1.375 


1.348 


1.392 


1.396 


E+ 






1.301 


1.256 


1.224 


1.277 


1.288 


E° 






1.364 


1.313 


1.288 


1.344 


1.345 




m 


4.45 ±0.22 


4.4581 


4.335 


4.757 


4.446 


4.650 


=0 






4.446 


4.315 


4.736 


4.440 


4.681 




[30] 


0.83 ±0.04 


0.8159 


0.7931 


0.8234 


0.8186 


0.8263 


=0 






0.8485 


0.8264 


0.8578 


0.8508 


0.8593 


!^ 


[22] 


0.62 ±0.09 


0.4499 


0.4906 


0.4414 


0.5165 


0.4299 


Q, 


[22] 


0.20 ±0.03 


0.1702 


0.1884 


0.1690 


0.1859 


0.1535 


A(1520) 


m 


1.57 ±0.44 


3.669 


2.724 


3.889 


3.382 


3.079 


e+(1540) 
9" (1540) 






2.68 
0.061 


2.41 
0.053 






5.71 
0.13 
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TABLE VI: Integrated multiplicities over limited rapidity windows around midrapidity obtained by fitting the measured 
distributions with single (G) or double (G+G) gaussians The yields at 10% and 20% centrality have been then multiplied 
by 1.08 and 1.32 respectively to convert them at 5% centrality trigger condition.. 



i Cli- LlVjlC 




R pfpfpTi r*p 




All — 9 


All — 1 


/Hnf 
X /oui 




^% 


yj 


G+G 


333.16 


176.82 


9.37 


K+ 


^% 




G+G 


57.16 


29.81 


3.33 


K" 


5% 


m 


G+G 


32.24 


16.90 


1.27 




5% 


m 


G 


4.327 


2.35 


0.04 


A 


10% 


[28j 


G 


22.56 


11.52 


0.49 


A 


10% 


[28j 


G 


3.039 


1.723 


0.58 




10% 




G 


2.75 


1.484 


1.2 




10% 


[30j 


G 


0.571 


0.3314 


0.71 


n 


20% 




G 


0.3095 


0.173 


0.73 


n 


20% 


[22J 


G 


0.126 


0.0789 


2.70 
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TABLE VII: Summary of fitted parameters in various heavy ion collisions at AGS and SPS in the framework of the SHM(7s) 
model. The 'common set' pararaeters have been obtained by fitting to the measured multiplicities of tt"*", K"*", K~, A, A and 
the participant nucleons in each collision. Also quoted minimum x^'s, the estimated radius of the EGG and the As parameter 
(sec Sect. 4). The re-scaled errors (see text) are quoted within brackets. For Pb-Pb at 30 A GeV of beam energy, we have 
constrained T and hb to lie on the fitted chemical freeze-out curve, as described in Sect. 4. 



Parameters 


Main analysis A 


Main analysis B 


Common set A 


Gommon set B 


Au-Au 11.6^ GeV 


T (MeV) 


118.1±3.5 (4.1) 


119.1±4.0 (5.4) 


119.2±2.1 (2.9) 


119.1±4.0 (6.6) 


AtB (MeV) 


555±12 (13) 


578±15 (21) 


556±12 (17) 


576.9±17.2 (29) 


7s 


0.652±0.069 (0.079) 


0.763±0.086 (0.12) 


0.645±0.042 (0.058) 


0.761±0.090 (0.15) 


VT^ exp[-0.7GeV/r] 


1.94±0.21 (0.24) 


1.487±0.18 (0.25) 


1.97±0.12 (0.17) 


1.494±0.21 (0.35) 


xVdof 


4.0/3 


5.5/3 


3.86/2 


5.5/2 


R (fm) 


9.31±0.69 (0.80) 


8.32±0.72 (0.97) 






As 


0.381±0.053 (0.061) 


0.490±0.084 (0.11) 


0.401±0.053 (0.074) 


0.487±0.089 (0.15) 


Pb-Pb 30A GeV 


T (MeV) 


139.5 


140.3 






/is (MeV) 


428.6 


428.7 






7s 


0.938±0.078 (0.13) 


1.051±0.103 (0.16) 






exp[-0.7GeV/T] 


6.03±0.50 (0.85) 


5.273±0.526 (0.80) 






xVdof 


5.75/2 


4.6/2 






As 


0.611±0.037 (0.062) 


0.683±0.086 (0.13) 






Pb-Pb AOA GeV 


T (MeV) 


147.6±2.1 (4.0) 


145.5±1.9 (3.5) 


148.6±2.1 (4.7) 


146.1±2.0 (4.0) 


HB (MeV) 


380.3±6.5 (13) 


375.4±6.4 (12) 


393±10 (22) 


390±10 (21) 


7s 


0.757±0.024 (0.046) 


0.807±0.025 (0.047) 


0.874±0.064 (0.14) 


0.961±0.079 (0.16) 


exp[-0.7GeV/T] 


8.99±0.37 (0.71) 


8.02±0.34 (0.63) 


8.09±0.55 (1.24) 


7.08±0.53 (1.1) 


xVdof 


14.7/4 


13.6/4 


10.1/2 


8.1/2 


R (fm) 


8.37±0.32 (0.61) 


8.37±0.31 (0.58) 






As 


0.507±0.025 (0.049) 


0.505±0.026 (0.048) 


0.621±0.064 (0.14) 


0.626±0.071 (0.14) 


Pb-Pb 80A GeV 


T (MeV) 


153.7±2.8 (4.7) 


151.9±3.4 (5.4) 


154.6±3.3 (7.2) 


152.2±3.5 (7.5) 


/is (MeV) 


297.7±5.9 (9.8) 


288.9±6.8 (11) 


300.7±9.4 (21) 


292.8±9.0 (19) 


7s 


0.730±0.021 (0.035) 


0.766±0.026 (0.042) 


0.741±0.057 (0.13) 


0.782±0.061 (0.13) 


VT^ exp[-0.7GeV/r] 


15.38±0.61 (1.0) 


14.12±0.65 (1.1) 


15.0±1.0 (2.3) 


13.7±0.95 (2.0) 


xVdof 


11.0/4 


10.4/4 


9.6/2 


9.3/2 


R (fm) 


9.03±0.41 (0.68) 


9.05±0.44 (0.71) 






As 


0.455±0.020 (0.034) 


0.461±0.020 (0.032) 


0.482±0.053 (0.12) 


0.4568±0.044 (0.095) 


Pb-Pb 158A GeV 


T (MeV) 


157.8±1.4 (1.9) 


154.8±1.4 (2.1) 


156.6±2.3 (3.3) 


152.7±2.1 (3.2) 


fiB (MeV) 


247.3±5.2 (7.2) 


244.5±5.0 (7.8) 


238.6±7.1 (10) 


232.4±7.7 (12) 


is 


0.843±0.024 (0.033) 


0.938±0.027 (0.042) 


0.722±0.053 (0.077) 


0.764±0.065 (0.097) 


VT^ exp[-0.7GeV/r] 


21.13±0.80 (1.1) 


18.46±0.69 (1.1) 


23.2±1.4 (2.0) 


21.1±1.4 (2.2) 


xVdof 


16.9/9 


21.6/9 


4.2/2 


4.5/2 


R (fm) 


9.41±0.26 (0.35) 


9.44±0.25 (0.39) 






A,,- 


0.500^0.018 (0.021) 


0.511^0.018 (0.028) 


0.120^0.037 (0.051) 


0.101^0.039 (0.058) 



21 



TABLE VIII: Fit results in Pb-Pb at 158^4 GeV with particle yields in limited rapidity window. 



Particle 


WA97 measured 


WA97 fitted 


NA49 measured 


NA49 fitted 


h" 


178±22 


157.2 






n~ 






176.8±9.8 


151.8 


K+ 






29.81±2.05 


30.63 


K" 






16.90±1.16 


18.24 


K°s 
</> 


21.9±2.4 


22.97 


2.35±0.34 


2.900 


A 


13.7±0.9 


13.75 


12.44±1.17 


15.35 


A 


1.8±0.2 


1.837 


1.86±0.13 


1.905 




1.5±0.1 


1.525 


1.603±0.079 


1.502 




0.37±0.06 


0.3782 


0.358±0.017 


0.3585 


n 






0.228±0.033 


0.1858 








0.104±0.016 


0.08894 




0.41±0.08 


0.3136 







TABLE IX: Fit results in Pb-Pb at 158^ GeV with particle yields in limited rapidity window. 



Parameters 


Ay=l 


Ay = 2 


WA97 central 


T (MeV) 
fXB (MeV) 

7s 

VT' (-xp[-().7Go\'/7T 


162.7±2.7 (5.1) 

229±12 (23) 
0.971±0.044 (0.083) 
5.55=0.31 (0.58) 


161.0±2.6 (4.2) 

223±13 (21) 
0.950±0.043 (0.070) 
iO. 71=0.59 (0.96) 


161.3±5.4 

218±19 
1.085±0.079 
1.73=0.52 




21.1/6 


16.0/6 


2.7/3 



TABLE X: Fit results in Pb-Pb at 158^ GeV with different models, as described in the text: SHM(SCV) (Strangeness Correla- 
tion Volume), SHM(TC) (Two Component model), SHM(75',7,) (light quaxk non-equilibrium model). Free fit parameters are 
quoted along with resulting minimum x^'s- The re-scaled errors (see text) are quoted withing brackets. For the SHM(7s,7q), 
the fit has been done by fixing 7^ = 1.6 near the absolute minimum (see Section 3.2). 



Parameters 


SHM(SCV) 


SHM(TC) 


SHM(7s,7,) 


T (MeV) 
HB (MeV) 
-fs 

VT^ exp[-0.7GeV/T] 
/ 


157.9±1.6 (3.3) 
261.5±2.4 (4.9) 

1.0 (fixed) 
18.62±0.52 (1.1) 
0.0253± 0.0067 (0.014) 


154.8±1.5 (1.9) 
237.1±7.0 (8.6) 
1.0 (fixed) 
15.50±0.54 (0.67) 

52.0±7.8 (9.6) 


140.4±1.1 (1.3) 
218.1±4.3 (5.2) 
0.929±0.027 (0.033) 
16.82±0.59 (0.72) 

1.6 (fixed) 


xVdof 


37.2/9 


13.7/9 


13.4/9 
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TABLE XI: Fitted parameters and multiplicities in pp collisions at a beam energy of 158 GeV, corresponding to ^/s = 17.2 
GeV. The re-scaled errors (see text) are quoted withing brackets. 



Parameter 


Value 


r (MeV) 




187.2±6.1 (9.3) 






5.79±0.85 (1.3) 


(ss) 


0.381±0.021 (0.032) 


xVdof 




16.1/7 




As 


0.224±0.019 (0.024) 


Particle 


Reference 


Measurement 


Fit 


7r+ 


[42] 


3.15±0.16 


3.257 


7r~ 


[42] 


2.45±0.12 


2.441 


tt" 






3.317 


K+ 


[42] 


0.21±0.02 


0.1901 


K~ 


[42] 


0.13±0.013 


0.09981 




[42] 


0.18±0.04 


0.1382 


V 






0.3918 


UJ 






0.3514 


t 


[29] 


0.012±0.0015 


0.01593 


V 






0.02576 


P+ 






0.4736 








0.3118 


P" 






0.4254 


K*+ 






0.07360 


K— 






0.02976 


K*° 






0.06192 


K*° 






0.03383 


P 






1.126 


P 


[42] 


0.040±0.007 


0.04364 


A++ 






0.2937 








0.007650 


A 


[43] 


0.115±0.012 


0.1123 


A 


[43] 


0.0148±0.0019 


0.01453 


E+ 






0.03480 


E- 






0.02310 


E° 






0.03004 


E" 






0.003317 


E+ 






0.004384 


E° 






0.003989 








0.001874 


1 — ^0 






0.002119 


=■+ 






0.0006902 








0.0006376 


n 






0.00003783 


n 






0.00002908 


A(1520) 


m 


0.012±0.003 


0.009155 


e+ 






0.005224 


6" 






0.0001515 
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FIG. 1: Above: measured versus fitted multiplicities in the statistical model supplemented with 7s parameter (SHM(7s)) in 
Au-Au collisions at a beam energy of 11.6^4 GeV in the fit A; also quoted the best-fit parameters. Below: residual distribution. 
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FIG. 2: Above: measured versus fitted multiplicities in the statistical model supplemented with 73 parameter in Pb-Pb collisions 
at a beam energy of 40A GeV in the fit A; also quoted the best-fit parameters. Below: residual distribution. 
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FIG. 3: Above: measured versus fitted multiplicities in the statistical model supplemented with 73 parameter in Pb-Pb collisions 
at a beam energy of 80A GeV in the fit A; also quoted the best-fit parameters. Below: residual distribution. 
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FIG. 4: Above: measured versus fitted multiplicities in the statistical model supplemented with 75 parameter in Pb-Pb collisions 
at a beam energy of 158A GeV in the fit A; also quoted the best-fit parameters. Below: residual distribution. Note that the 
A(1520) was not used in the fit (see text). 
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FIG. 5: Comparison between measured and calculated (in fits A and B) (K"'")/(7r"'")and (7r~)/(Np)ratios as a function of the 

centre-of-mass energy in the examined collisions. For the SPS energy points the statistical errors are indicated with solid lines, 
while the contribution of the common systematic error is shown as a dotted line The lines connect the fitted values. 



Oh 



2.00 
1.80 
1.60 
1.40 
1.20 
1.00 
0.80 
0.60 
0.40 
0.20 
0.30 
0.28 
0.26 
0.24 
0.22 
0.20 
0.18 
0.16 
0.14 
0.12 
0.10 







exp. <7i:'>/<Np> 
th. <7i">/<Np> A 
th. <7i">/<Np> B 



exp. <K^>/<7i;'> 
th. <K+>/<7i^> A 
th. <K"'>/<7t''> B 




10 



15 
V S 



20 



NN 



(GeV) 



28 



FIG. 6: Spectrum of known light-flavoured hadronic species up to a mass of 1.8 GeV. 
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FIG. 7: Left: primary yields of various particles as a function of the cut-off on the hadronic mass spectrum. Right: fitted 7s, 
baryon-chemical potential and temperature as a function of the cut-off on the hadronic mass spectrum. 
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FIG. 9: Above: measured versus fitted multiplicities in the statistical model supplemented with strangeness suppression in pp 
collisions at a beam energy of 158 GeV corresponding to \/i = 17.2 GeV; also quoted the best-fit parameters. Below: residual 
distribution. 
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FIG. 10: Minimum of multiplicity fits in Pb-Pb collisions at a beam energy of 158^ GeV as a function of a fixed light quark 

non equilibrium parameter 7,. The round dots are the x^'s obtained with the main data sample, whilst triangular dots are 
those obtained with pion multiplicities lowered by one standard deviation and all others unchanged. The vertical dashed line 
indicates the condensation point. 
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FIG. 11: Chemical freeze-out points in the /is — T plane in various heavy ion collisions. The full round dots refer to Au-Au 
at 11.6 and Pb-Pb collisions at 40, 80, 158A GeV obtained in the analysis A, whilst the hollow square dot has been obtained 
in ref. fTH by using particle ratios measured at midrapidity in Au-Au collisions at %/sjvjv ~ ^"^0 GeV. The hollow round dot 
without error bars refers to Pb-Pb collisions at 30^ GeV and has been obtained by forcing T and ^lb to lie on the parabola 
fitted to the full round dots. 
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FIG. 12: Strangeness non-equilibrium parameter 7s as a function of the nucleon-nucleon centre-of-mass energy. Full dots refer 
to fit A, hollow dots to fit B. 



1.2 - 



0.8 - 




0.6 - 



0.4 - 



0.2 - 



2 4 6 8 10 12 14 16 18 20 



V s.... (GeV) 



35 



FIG. 13: Measured {K+)/(7r+)ratio as a function of the fitted baryon-chemical potential. The full square dot is a preliminary 
full phase space measurement in Au-Au collisions at \/sjvjv ~ 200 GeV i23] a^nd the error is only statistical; the arrow on the 
left signifies that its associated baryon chemical potential is lower than that estimated at \/sjvjv ~ ^'^^ GeV used here. 
For the SPS energy points the statistical errors are indicated with solid lines, while the contribution of the common systematic 
error is shown as a dotted line. Also shown the theoretical values for a hadron gas along the fitted chemical freeze-out curve 
shown in fig. 1111 for different values of 7s. 
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FIG. 14: As estimated from the fits A (full dots) and B (hollow dots) as a function of the fitted baryon-chemical potential. Also 
shown the theoretical values for a hadron gas along the fitted chemical freeze-out curve shown in fig. 1111 for different values of 
7s- 
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